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ABSTRACT 


The  measurement  of  regional  attenuation  QA  can  produce  method  dependent  results.  The  discrepancies  among 
methods  are  due  to  differing  parameterizations  (e.g.,  geometrical  spreading  rates),  employed  data  sets  (e.g.,  choice 
of  path  lengths  and  sources),  and  methodologies  themselves  (e.g.,  measurement  in  the  frequency  or  time  domain). 
We  apply  the  coda  normalization  (CN),  two-station  (TS),  reverse  two-station  (RTS),  source-pair/receiver-pair 
(SPRP),  and  the  new  coda-source  normalization  (CS)  methods  to  measure  Q  of  the  regional  phase,  Lg  (Qi.o).  and  its 
power-law  dependence  on  frequency  of  the  form  Qtf*  with  controlled  parameterization  in  the  well-studied  region  of 
northern  California  using  a  high-quality  data  set  from  the  Berkeley  Digital  Seismic  Network.  We  test  the  sensitivity 
of  each  method  to  changes  in  geometrical  spreading,  Lg  frequency  bandwidth,  the  distance  range  of  data,  and  the  Lg 
measurement  window.  For  a  given  method,  there  are  significant  differences  in  the  power-law  parameters,  Qo  and  //. 
due  to  perturbations  in  the  parameterization  when  evaluated  using  a  conservative  pairwise  comparison.  The  CN 
method  is  affected  most  by  changes  in  the  distance  range,  which  is  most  probably  due  to  its  fixed  coda  measurement 
window.  Since,  the  CS  method  is  best  used  to  calculate  the  total  path  attenuation,  it  is  very  sensitive  to  the 
geometrical  spreading  assumption.  The  TS  method  is  most  sensitive  to  the  frequency  bandwidth,  which  may  be  due 
to  its  incomplete  extraction  of  the  site  term.  The  RTS  method  is  insensitive  to  parameterization  choice,  whereas  the 
SPRP  method  as  implemented  here  in  the  time-domain  for  a  single  path  has  great  error  in  the  power-law  model 
parameters  and  77  is  greatly  affected  by  changes  in  the  method  parameterization.  When  presenting  results  for  a  given 
method  it  is  best  to  calculate  Qqf  for  multiple  parameterizations  using  some  a  priori  distribution.  We  also 
investigate  the  difference  in  power-law  Q  calculated  among  the  methods  by  considering  only  an  approximately 
homogeneous  subset  of  our  data.  All  methods  return  similar  power-law  parameters,  though  the  95%  confidence 
region  is  large.  We  adapt  the  CS  method  to  calculate  (9/„  tomography  in  northern  California.  Preliminary  results 
show  that  by  correcting  for  the  source,  tomography  with  the  CS  method  may  produce  better  resolved  attenuation 
structure. 
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OBJECTIVES 


Understanding  of  regional  attenuation  Q1  can  help  with  structure  and  tectonic  interpretation  (Aleqabi  and 
Wysession,  2006;  Benz  et  ah,  1997;  Frankel,  1990),  and  correcting  for  the  effects  of  attenuation  can  lead  to  better 
discrimination  of  small  nuclear  tests  (e.g.  Baker  et  ah,  2004;  Mayeda  et  ah,  2003;  Taylor  et  ah,  2002).  Present 
threshold  algorithms  for  event  identification  rely  on  Q  models  that  are  derived  differently,  and  the  models  can  vary 
greatly  for  the  same  region.  For  example,  recent  one-dimensional  (l-D)  Q  studies  in  South  Korea  find 
frequency-dependent  QLg  that  at  1  FIz  range  from  450  to  900  (Chung  and  Lee  2003;  Chung  et  ah,  2005).  Another 
example  is  the  case  of  Tibet  where  a  wide  variety  of  Q  values  have  been  reported  (e.g.  Fan  and  Lay,  2003;  Xie  et  ah 
2004).  In  order  to  reliably  use  reported  Q  estimates  for  either  monitoring  applications,  or  for  tectonic  interpretation  it 
is  essential  to  know  the  uncertainty  in  the  estimate.  Commonly,  individual  studies  will  present  aleatoric  (random) 
uncertainty,  however  epistemic  (bias)  uncertainty  is  not  possible  to  assess,  when  only  a  single  method  and 
parameterization  is  considered.  In  order  to  better  understand  the  effects  of  different  methods  and  parameterizations 
on  Q  models,  we  implement  four  popular  methods  and  one  new  method  to  measure  Q  of  the  regional  seismic  phase, 
Lg  (QLg),  using  a  high-quality  data  set  from  the  Berkeley  Digital  Seismic  Network  (BDSN).  The  CN  method  is 
implemented  in  the  time  domain  for  paths  leading  to  a  common  station  and  it  returns  a  stable  Q  measurement  when 
the  region  near  a  station  is  homogenous.  The  CS  method  uses  previously  calculated  coda-derived  source  spectra  to 
remove  the  source  term  in  the  frequency  domain  and  is  best  suited  to  calculate  an  effective  Q  for  a  given  path.  The 
TS  and  RTS  methods  are  implemented  in  the  frequency  domain  and  the  calculated  Q  is  more  stable  due  to  the 
extraction  of  the  source  term.  The  RTS  method  produces  a  power-law  Q  with  less  error  than  the  TS  method  due  to 
its  additional  extraction  of  the  site  terms,  though  it  is  more  restrictive  in  its  data  requirements.  The  SPRP  method  is 
the  RTS  method  with  a  relaxation  of  the  data  requirements  and  is  implemented  in  the  time  domain  here. 

Through  this  approach  we  identify  both  aleatoric  and  epistemic  uncertainty.  With  a  more  complete  knowledge  of 
uncertainty  it  will  be  possible  to  better  assess  the  results  of  published  attenuation  studies  and  the  presented 
multi-method  analysis  procedure  employed  in  future  efforts  can  lead  to  improved  estimates  of  regional  Q. 

RESEARCH  ACCOMPLISHED 


The  data  set  consists  of  158  earthquakes  recorded  at  16 
broadband  (20  sps)  three-component  stations  of  the 
BDSN  between  1992  and  2004  (Figure  1).  The  wide 
distribution  of  data  parameters  allows  for  sensitivity 
testing.  We  calculate  QLg  by  fitting  the  power-law 
model,  Qof'  using  five  different  methods.  The  first  two 
methods  use  the  seismic  coda  to  correct  for  the  source 
effect.  The  last  three  methods  use  a  spectral  ratio 
technique  to  correct  for  source,  and  possibly  site 
effects.  In  the  following  we  summarize  the  methods 
and  point  out  significant  differences.  Our  philosophy  in 
presenting  each  of  the  methods  is  to  maintain  the 
approach  and  style  of  the  popular  version  of  each 
method  as  close  as  possible.  Later,  we  will  attempt  to 
normalize  each  of  the  methods  for  comparison  and 
sensitivity  testing. 

Coda  Normalization  (CN) 

The  CN  method  uses  the  local  shear-wave  coda  as  a 
proxy  for  the  source  and  site  effects,  thus  amplitude 
ratios  remove  these  two  effects  from  the  5-wave 
spectrum  (Aki,  1980;  Yoshimoto  et  ah,  1993).  In  his 
236  238  240  242  original  application,  Aki  (1980)  assumed  that  the  local 

shear-wave  coda  was  homogeneously  distributed  in 

Figure  1.  Events  (stars),  stations  (inverted  triangles),  space  and  time.  For  the  current  study  region.  Figure  1 
used  in  the  study.  of  Mayeda  et  al.  (2005)  shows  that  the  coda  at  ~1  Hz  is 
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in  fact  homogeneous,  at  least  up  to  ~240  km.  This  method  assumes  the  Lg  amplitude  ALg  at  a  given  distance  r  and 
frequency/ can  be  estimated  by 


AAf,r)  =  S  {f)R  {9)1  (f)P(f)G(r)  exp 


-rtf 

QU 


(1) 


where  S(f)  is  the  source  spectrum  and  R(  6)  is  the  source  radiation  in  the  source-receiver  direction  6.  P(J)  is  the  site 
term,  /(/)  is  the  instalment  term,  and  G(r)  is  the  geometrical  spreading  term,  approximated  here  as  an  inverse 
power-law  where  yis  the  spreading  rate  and  is  given  in  Table  1.  The  final  term  is  an  apparent  attenuation,  where  U 
is  the  Lg  group  velocity,  which  is  fixed  at  3.5  km/s  for  this  and  all  other  methods.  The  CN  method  also  assumes  that 
the  coda  spectrum  C(f)  is  approximately  equal  to  the  source  spectrum  at  a  given  critical  propagation  time  tc-  The 
coda  excitation  term  is  assumed  to  be  constant  at  all  distances  for  a  given  tc.  If  the  source  radiation  is  smoothed  by 
considering  several  sources  at  many  source-receiver  directions  we  can  take  the  ratio  of  ALg  to  C,  measured  at  tc, 
which  effectively  removes  instillment,  site,  and  source  contributions  resulting  in  only  the  geometrical  spreading  and 
attenuation  terms.  The  natural  log  of  this  spectral  ratio  taken  at  discrete  frequency  bands  (between  0.25,  0.5,  1,  2,  4, 
and  8  Hz)  results  in  the  equation  of  a  line  as  a  function  of  distance  and  the  slope  is  related  to  Q1 .  Q'x  at  the  center 
frequency  of  each  band  then  reveals  a  power-law  model  for  each  station. 

ALg  is  the  maximum  envelope  amplitude  in  each  bandpassed  (8-pole  acausal  Butterworth  filter),  windowed 
(according  to  the  window  parameter  in  Table  1)  and  tapered  raw  vertical  trace.  C  is  the  root-mean-square  (rms) 
amplitude  in  each  bandpassed  10  second  window  centered  on  a  tc  of  150  sec.  Data  were  excluded  if  either  ALg  or  C 
had  a  SNR  less  than  two,  where  noise  is  measured  as  the  maximum  amplitude  in  a  window  the  same  length  as  ALg 
prior  to  the  event.  This  method  is  similar  to  that  of  Chung  and  Lee  (2003),  whereas  Frankel  (1990)  used  a  weighted 
average  of  the  smoothed  coda  to  measure  C.  We  calculate  Q'x  with  all  records  at  a  given  station,  where  the  slope  is 
calculated  with  an  iteratively  weighted  least-squares  method.  The  resulting  QA  are  then  fit  in  the  log  domain  as  a 
function  of  midpoint  frequency  with  a  weighted  (the  squared  inverse  of  the  standard  error  in  each  Q  ]  measurement) 
least-squares  line  to  calculate  the  power-law  parameters. 

Coda-Source  Normalization  (CS) 

The  CS  method  uses  the  stable,  coda-derived  source  spectra  to  isolate  the  path  attenuation  component  of  the  Lg 
spectmm  (Walter  et  ah,  2007).  This  method  assumes  ALg  is  represented  as  in  equation  (1)  with  S(f)  described  as  in 
Aki  and  Richards  (2002),  G(r)  is  a  critical  distance  formulation  (Street  et  ah,  1975).  We  assume  a  site  term  P(f)  of 
unity  and  thus  any  site  effect  is  projected  into  the  path  attenuation  term. 

The  windowed  (according  to  the  window  parameter  in  Table  1)  and  tapered  transverse  component  is  transferred  to 
velocity  and  its  Fourier  amplitude  is  calculated.  ALg  is  then  the  mean  of  the  Fourier  amplitude  for  fixed  discrete 
frequency  bands  (between  0.2,  0.3,  0.5,  0.7,  1,  1.5,  2,  3,  4,  6,  and  8  Hz).  Path  attenuation  can  then  be  extracted  with 
the  log  transform  where  the  same  frequency  bands  are  used  to  calculate  the  source  spectra,  S(f),  and  P(J)  is  fixed  to 
unity.  Source  spectra  derived  from  the  coda  are  calculated  via  the  methodology  of  Mayeda  et  ah  (2003)  and  from  the 
northern  California  study  of  Mayeda  et  ah  (2005).  Q(j)  is  only  calculated  for  records  where  ALg  is  two  times  the 
amplitude  of  the  pre-event  signal  (SNR  >  2).  Q  at  the  center  frequency  of  each  band  then  reveals  a  power-law  model 
for  each  event-station  path.  We  fit  a  least-squares  line  in  the  log  domain  (a  robust  regression  gave  similar  results) 
and  the  intercept  term  is  then  the  log  transform  of  Q0  and  the  slope  is  //. 


Table  1.  Method  parameterization. 


Group 

Spreading  exponent 

[y] 

Measurement  band 
(Hz) 

Epicentral  distance 
M  (km) 

Lg  Velocity  window 
(km/s) 

Control 

0.50 

0.50-8 

100  -400 

2. 6-3. 5 

Test  1  (y) 

Test  2  (Bandwidth) 
Test  3  (Distance) 

Test  4  (Window) 

0.83 

0.25-4 

100-700 

3. 0-3. 6 
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Two-Station  (TS) 

The  TS  method  takes  the  ratio  of  Lg  recorded  at  two 
different  stations  along  the  same  narrow  path  from  the 
same  event  in  order  to  remove  the  common  source  term 
(e.g.,  Chavez  and  Priestley,  1986;  Xie  and  Mitchell, 
1990).  We  implement  this  method  in  the  frequency 
domain  and  take  the  ratio  of  two  terms  with  the  form  of 
equation  (1),  which  can  then  be  transformed  to  the 
log-domain  and  a  linear  regression  is  possible  to 
calculate  the  power-law  parameters.  However,  random 
error  due  to  propagation  can  produce  a  negative 
argument  in  the  exponential  term  of  equation  (1)  at 
some  frequencies  (Xie,  1998),  which  prohibits  analysis 
in  the  log-domain.  Therefore,  we  perform  a  non-linear 
regression  on  the  argument  that  minimizes  the  sum  of 
squares  error  on  the  power-law  function  in  the 
least-squares  sense  (Bates  and  Watts,  1988). 

Reverse  Two-Station  (RTS) 

The  RTS  method  uses  two  TS  setups,  where  a  source  is 
on  either  side  of  the  station  pair  in  a  narrow  azimuthal 
window  (Chun  et  al.,  1987).  The  two  ratios  are 
combined  to  remove  the  common  source  and  site  terms. 

Source-Pair/Receiver-Pair  (SPRP) 

The  SPRP  method  is  the  RTS  method  with  a  relaxation 
on  the  narrow  azimuthal  window  requirement  (Shih 
et  al.,  1994).  We  implement  this  method  in  the 
time-domain.  Unlike  the  RTS  method,  data  are  no 
longer  restricted  by  a  given  azimuth  but  by  a  distance 
formulation.  ALg  is  the  maximum  zero-to-peak 
amplitude  in  each  bandpassed  (8-pole  acausal 
Butterworth  filter),  windowed  (according  to  the  window 
parameter  in  Table  1)  and  tapered  vertical  component 
record  that  has  been  transferred  to  velocity.  The 
equation  is  least-squares  fit  as  a  function  of  the  effective 
interstation  distance  for  the  same  discrete  frequency 
bands  as  in  the  CN  method,  where/is  the  midpoint  of 
these  frequency  bands.  The  slope  of  the  fit  is  a  function 
of  (9  1  in  the  band  that  was  measured.  The  resulting  Q'1 
are  then  fit  in  the  log  domain  as  a  function  of  midpoint 
frequency  with  a  weighted  (the  squared  inverse  of  the 
standard  error  in  each  QA  measurement)  least-squares 
line  to  calculate  the  power-law  parameters. 

Method  Comparison 

Since  each  method  has  a  different  data  requirement  it  is 
improper  to  compare  the  methods  with  the  full  data  set. 
For  example,  the  CN  method  will  sample  geology  at  all 
back-azimuths  relative  to  a  station,  whereas  the  RTS 
method  is  restricted  to  a  narrow  azimuthal  window 
aligned  roughly  along  a  pair  of  stations  and  events. 


Figure  2.  Method  comparison,  (a)  Map  (same  region 
as  Figure  1)  of  the  subset  used  in  the 
comparison  analysis.  Data  are  in  a  small 
region  near  the  San  Francisco  Bay  Area, 
primarily  along  the  Franciscan  block. 

(b)  Power-law  parameters  and  their 
empirical  95%  confidence  regions  are 
given.  The  intersecting  region  is  shaded 
grey. 
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QO 

Figure  3.  Method  comparison  with  tomographic 
results  of  (Mayeda  et  al.  2005).  (a)  Region 
where  2-D  direct  wave  attenuation 
coefficients  are  used,  which  covers  the  same 
area  as  the  paths  and  stations  in  Figure  2a. 
(b)  Comparison  of  tomographic  results 
where  empirical  distribution  (light  grey)  is 
from  data  at  each  node  with  1-D  results. 


In  an  attempt  to  normalize  the  data  set  used  for  each 
method,  we  restrict  the  data  to  lie  in  a  small  region 
along  the  Franciscan  block  (Figure  2a).  We  implement 
all  five  methods  to  calculate  Qtf1  in  the  region 
(Figure  2b).  The  populations  are  then  smoothed  with  a 
two-dimensional  gaussian  kernel  (Venables  and  Ripley, 
2002)  to  produce  an  empirical  distribution  so  that  the 
95%  confidence  region  can  be  estimated.  The  grey 
region  in  Figure  2b  represents  a  parameter  space  that 
fits  all  studies.  Mayeda  et  al.  (2005)  present  Q 
tomography  for  Northern  California  in  order  to  compare 
1-D  and  2-D  methods  to  calculate  both  coda  and  direct 
wave  (5,  Lg,  or  surface  wave)  attenuation.  We  extend 
the  analysis  for  comparison  with  the  results  from  the 

1- D  analysis  of  the  sub-region.  Power  law  parameters 
from  the  Mayeda  et  al.  (2005)  study  are  calculated  by 
fitting  a  least-squares  line  to  the  Q  estimated  for  each 
frequency  band  at  the  midpoint  of  the  band  in  the  log 
domain.  We  extract  the  power-law  parameters  at  points 
within  the  sub-region  (Figure  3a)  and,  as  above,  we 
produce  an  empirical  distribution  (Figure  3b).  The  range 
in  rj  and  variance  of  Q0  are  similar  between  the  1-D  and 

2- D  results,  but  the  mean  of  the  Q0  distribution  is 
shifted  by  about  30.  This  may  be  due  to  some 
regularization  effects.  This  analysis  shows  that  some  of 
the  variability  in  the  1-D  analysis  is  due  to  2-D 
structure. 

Sensitivity  Tests 

Using  the  complete  data  set,  we  investigated  how  the 
choice  of  parameterization  affects  the  results.  In  each 
test,  only  one  parameter  was  varied,  and  Qtf'  was 
calculated  with  each  of  the  methods.  The  varied 
parameters  are  geometrical  spreading  rate,  measurement 
bandwidth,  epicentral  distance,  and  the  Lg  window.  The 
values  of  the  varied  parameters  are  listed  in  Table  1, 
where  the  range  was  chosen  based  on  the  values  used  in 
previous  studies. 

For  the  CN  method,  standard  error  regions  were 
constructed  from  the  covariance  of  the  power-law 
model  parameters  estimated  by  bootstrapping  the 
residuals  of  the  weighted  least-squares  fit  1000  times 
(Aster  et  al.,  1996).  Figure  4a  shows  the  standard  error 
regions  for  each  test  at  station  PKD.  All  tests  cluster 
around  the  control  parameters  except  the  distance  test 
(Test  3).  To  assess  the  significance  of  model 
parameterization  differences  we  perform  an  analysis  of 
covariance  (ANCOVA)  for  the  weighted  least-squares 
regression  with  Tukey’s  honest  significant  difference 
pairwise  comparison  tests  (Faraway,  2004).  A 
difference  in  the  model  parameters  is  only  significant  if 
the  95%  confidence  region  of  the  mean  difference  in  the 
model  parameters  between  two  tests  does  not  include 
zero.  We  group  all  significant  differences  between  a 
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given  test  and  the  control  parameterization  and  plot  the 
median,  25th  and  75th  percentile  values  of  that  group 
(Figure  4b).  In  this  way,  we  can  try  and  separate  aleatoric 
uncertainty  due  to  poorly  constrained  power-law  model 
parameters  and  epistemic  uncertainty  due  to  the  choice  of 
parameterization  for  each  method,  and  one  can  think  of  the 
confidence  regions  in  panel  a)  of  Figures  4  as  the  aleatoric 
uncertainty,  and  the  values  in  panel  b)  as  epistemic 
uncertainty.  There  is  a  significant  difference  for  almost  all 
CN  method  comparisons  in  rj,  and  the  greatest  difference  for 
both  model  parameters  is  when  the  epicentral  distance  of  the 
data  set  is  changed  (Test  3).  This  is  due  to  the  fixed  time  tc 
at  which  the  coda  is  measured,  where  for  greater  distances  it 
is  more  appropriate  to  increase  to . 

Standard  error  regions  and  pairwise  comparisons  are 
calculated  for  the  CS  method  as  described  above,  though  the 
residuals  and  ANCOVA  are  for  a  direct  linear  regression. 
For  most  tests  only  a  small  fraction  of  the  comparisons  are 
significant.  However,  when  y  is  changed  in  Test  1,  there  is  a 
significant  difference  in  Q0  for  39%  of  the  path 
comparisons,  where  the  median  difference  is  almost  50.  This 
effect  highlights  the  difficulty  in  extracting  an  intrinsic  Q 
from  the  full  path  attenuation  when  examining  a  single  path. 
The  CS  method  is  best  for  evaluating  the  total  path  term. 

Since  the  TS  and  RTS  methods  require  nonlinear 
regressions,  we  estimate  covariance  matrices  from  the 
bootstrapped  power-law  model  parameter  populations. 
ANCOVA  is  performed  with  this  estimated  covariance  and 
the  pairwise  comparisons  are  made  with  the  results.  A 
change  in  epicentral  distance  does  not  significantly  affect 
the  power-law  parameters  for  both  the  TS  and  RTS 
methods,  but  a  change  in  bandwidth  (Test  2)  produces  an 
interquartile  range  of  0.05  to  0.22  for  the  difference  in  // 
using  the  TS  method.  The  TS  method  is  sensitive  to  site 
effects  and  this  difference  may  be  due  to  site  effects  that  are 
different  below  1  Hz  than  they  are  above  it.  For  several 
stations  in  the  BDSN  this  seems  to  be  the  case  (Malagnini 
et  al.,  2007).  The  RTS  method  doesn’t  suffer  from  this  same 
dependency  and  its  median  significant  differences  are  low 
for  all  tests. 

As  previously  stated,  the  SPRP  method  implemented  in  the 
time  domain  requires  a  distribution  of  effective  interstation 
distances  that  can  best  be  given  when  several  interstation 
paths  are  considered.  However,  it  should  be  able  to 
constrain  Q^f1  for  a  single  interstation  path,  and  in  order  to 
allow  for  comparison  with  the  implementation  of  the  other 
interstation  methods,  TS  and  RTS,  we  carry  out  the  method 
on  an  interstation  basis.  Due  to  such  large  standard  error 
regions  only  around  half  of  the  pairwise  comparisons  give  a 
significant  difference  in  Q0.  However,  the  same 
comparisons  reveal  a  large  difference  in  ij  for  all  but  the  y 
test  (Test  1). 
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Figure  4.  Parameterization  effects  of  the  coda- 
normalization  method,  a)  Power-law 
parameters  (Q0,  rj)  for  each  choice  of 
parameterization  and  the  standard  error 
region,  b)  Results  of  significant 
difference  in  pairwise  comparisons 
between  the  Control  parameterization 
and  its  deviations  (similar  symbol  as  a) 
for  all  measurements  in  the  method.  The 
upper  right  box  gives  percentage  of 
measurements  that  had  a  significant 
difference  and  the  symbols  are  at  the 
median  difference  (SQ0,  8rj)  with  upper 
(3rd  quartile)  and  lower  (1st  quartile) 
bounds  given  by  the  bars. 
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Figure  5.  Attenuation  tomography  using  amplitudes  (a)  where  the  source  and  site  term  is  solved  for  in  the 
inversion  and  (b)  where  the  source  term  is  removed  from  the  amplitudes  using  a  coda-derived 
moment  rate. 


Discussion 

Each  method  analyzed  here  is  employed  for  different  types  of  investigation.  Each  method  has  different  advantages, 
disadvantages  and  assumptions.  The  CN  method  returns  a  stable  Q  measurement  when  the  region  near  a  station  is 
homogenous.  The  CS  method  is  best  suited  to  calculate  an  effective  Q  for  a  given  path,  where  the  site  term  is 
mapped  into  the  path  attenuation.  Also,  since  it  measures  the  path  directly  from  the  event  to  station,  there  is  a 
trade-off  between  geometrical  spreading  and  effective  Q.  If  the  uncertainties  in  the  type  of  geometrical  spreading  are 
large,  then  it  may  be  best  to  test  several  forms  of  spreading,  or  to  fold  the  spreading  term  into  the  entire  path  effect  if 
this  is  appropriate  for  the  application.  The  CS  method  can  be  used  to  calculate  corrected  amplitudes  for  use  in  a 
tomographic  inversion.  We  have  created  such  a  scheme  by  adapting  the  method  of  Phillips  et  al.  (2005)  and 
preliminary  results  show  that  this  method  may  resolve  structure  more  tightly  (Figure  5). 

The  TS  and  RTS  methods  are  more  stable  due  to  the  extraction  of  the  source  term.  The  RTS  method  produces  the 
least  error  due  to  its  additional  extraction  of  the  site  terms,  though  it  is  more  restrictive  in  its  data  requirements.  Xie 
(2002)  calculates  the  bias  due  to  the  site  term  assumption  in  the  TS  method  and  finds  that  it  is  small.  In  order  to  test 
this  assumption  we  compare  the  average  power-law  parameters  for  paths  calculated  by  both  the  TS  and  RTS 
methods  (Figure  6a).  The  values  of  the  parameters  are  approximately  the  same  for  both  methods,  though  there  is 
scatter.  A  more  direct  test  is  to  compare  the  power-law  parameters  calculated  for  paths  to  station  BKS  and  new  data 
from  a  nearly  co-located  BRK  (Figure  6b).  Malagnini  et  al.  (2007)  find  a  significant  difference  in  the  site  term 
between  BKS  and  BRK  and  this  difference  is  evident  in  Figure  6b.  Stacking  ratios  with  common  interstation  paths 
could  reduce  the  variance,  but  this  is  only  appropriate  for  tectonically  stable  areas.  Aster  et  al.  (1996)  calculates 
spectral  ratios  with  the  multi-taper  method  and  is  able  to  produce  more  stable  spectra  and  a  more  realistic  variance  in 
the  spectral  measurement. 
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The  SPRP  method  is  the  RTS  method  with  a  relaxation 
of  the  data  requirements  and  is  appropriate  for  very 
laterally  homogeneous  Q.  The  SPRP  method  is 
implemented  in  the  frequency  domain  by  Fan  and  Lay 
(2003),  and  in  the  time  domain  by  Shih  et  al.  (1994)  and 
Chung  et  al.  (2005),  where  they  find  clusters  in  small 
regions  that  are  very  different  from  the  overall  1-D  Q 
model.  The  SPRP  method  in  the  time  domain  is  much 
better  suited  for  a  large  homogeneous  region,  where 
several  interstation  regions  can  be  grouped  together.  In 
the  implementation  here,  we  calculate  Qi/]  for  each 
interstation  path  that  fits  the  above  criteria  (<41%  of  the 
available  paths),  which  results  in  pooling  of  data  points 
near  the  true  interstation  distance.  This  can  greatly 
effect  the  linear  regression  and  produce  large  error  in 
the  model  parameters. 

CONCLUSIONS  AND  RECOMMENDATIONS 

We  apply  the  coda  normalization  (CN),  two-station 
(TS),  reverse  two-station  (RTS),  source-pair/ 
receiver-pair  (SPRP),  and  the  new  coda-source 
normalization  (CS)  methods  to  measure  QLg  and  its 
power-law  dependence  (Qt/1)  in  northern  California  in 
order  to  understand  the  variability  due  to 
parameterization  choice  and  method.  We  investigate  the 
reliability  of  the  methods  by  comparing  them  with  each 
other  for  an  approximately  homogeneous  region  in  the 
Franciscan  block  near  the  San  Francisco  Bay  Area.  All 
methods  return  similar  power-law  parameters, 
especially  in  their  95%  confidence  regions.  If  we 
consider  the  joint  distributions  of  each  method  Q0  =  85 
±  40  and  ;/=  0.65  ±  0.35  (both  -95%  Cl),  where  ^  is 
not  as  well  constrained.  We  test  the  sensitivity  of  each 
method  to  changes  in  geometrical  spreading,  Lg 
frequency  bandwidth,  the  distance  range  of  data,  and  the 
Lg  measurement  window.  For  a  given  method,  there  are 
significant  differences  in  the  power-law  parameters,  Qo 
and  //.  due  to  perturbations  in  the  parameterization  when 
evaluated  using  a  conservative  pairwise  comparison. 

The  CN  method  is  affected  most  by  changes  in  the 
distance  range,  which  is  most  probably  due  to  its  fixed 
coda  measurement  window  or  the  fact  that  at  larger 
distances  the  coda  is  not  homogeneously  distributed. 

Since,  the  CS  method  is  best  used  to  calculate  the  total 
path  attenuation,  it  is  very  sensitive  to  the  geometrical 
spreading  assumption.  The  TS  method  is  most  sensitive 

to  the  frequency  bandwidth,  which  may  be  due  to  its  incomplete  extraction  of  the  site  term.  The  RTS  method  is 
insensitive  to  parameterization  choice,  whereas  the  SPRP  method  as  implemented  here  in  the  time -domain  for  a 
single  path  has  great  error  in  the  power-law  model  parameters  and  //  is  greatly  affected  by  changes  in  the  method 
parameterization.  When  presenting  results  for  a  given  method  it  is  best  to  calculate  Qi/]  for  multiple 
parameterizations  using  some  a  priori  distribution.  We  plan  to  implement  the  methods  introduced  here  to  find  1-D 
and  2-D  attenuation  in  the  Y ellow  Sea  /  Korean  Peninsula  region. 


Figure  6.  An  investigation  into  site  effects. 

(a)  Power-law  parameters  for  paths 
measured  by  both  the  RTS  and  TS 
methods,  (b)  Power-law  parameters 
measured  at  nearly  co-located  stations 
BKS  and  BRK  using  the  TS  method. 
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